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Abstract 

We consider the hypothesis that Dark Matter (DM) has tree-level interactions only with leptons. 
Such a framework, where DM recoils against electrons bound in atoms, has been proposed as an 
explanation for the annually modulated scintillation signal in DAMA/LIBRA versus the absence 
of a signal for nuclear recoils in experiments like CDMS or XENON10. However, even in such a 
leptophilic DM scenario there are loop induced DM-hadron interactions, where photons emitted 
from virtual leptons couple to the charge of a nucleus. Using a general effective field theory 
approach we show that, if such an interaction is induced at one or two-loop level, then DM-nucleus 
scattering dominates over DM-electron scattering. This is because the latter is suppressed by the 
bound state wave function. One obtains a situation similar to standard DM-nucleus scattering 
analyses with considerable tension between the results of DAMA and CDMS/XENON10. This 
conclusion does not apply in the case of pseudoscalar or axial vector coupling between DM and 
leptons, where the loop diagrams vanish. In this case the explanation of the DAMA signal in terms 
of DM-electron scattering is strongly disfavored by the spectral shape of the signal. Furthermore, if 
DM can annihilate into neutrinos or tau leptons, the required cross sections are excluded by many 
orders of magnitude using the Super-Kamiokande bound on neutrinos from DM annihilations in 
the Sun. 
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I. INTRODUCTION 



The DAMA collaboration provided strong evidence for an annually modulated signal in 
the scintillation light from sodium iodine detectors. The combined data from DAMA/Nal H 
(7 annual cycles) and DAMA/LIBRA p[ (4 annual cycles) with a total exposure of 0.82 ton yr 
show a modulation signal with 8.2a significance. The phase of this modulation agrees 
with the assumption that the signal is due to the scattering of Weakly Interacting Massive 
Particles (WIMPs) forming the Dark Matter (DM) halo of our Galaxy. 

However, many interpretations of this signal in terms of DM scattering are in conflict 
with constraints from other DM direct detection experiments. Spin-independent elastic 
WIMP-nucleon scattering accounting for the DAMA modulation is tightly constrained by 
bounds from several experiments, most notably from CDMS J3] and XENON10 M. While 
conventional WIMPs with masses m x > 50 GeV are excluded by many orders of magnitude, 
light WIMPs with < 10 GeV masses might be marginally compatible with the constraints, 
see, e.g. |7|, |S], |S], |TD], [H], [T2j for recent studies. Spin-dependent couplings to protons 

can account for the DAMA signal without being in conflict with CDMS and XENON10, 
but in this case strong constraints from COUPP JlJ), KIMS |L4|, PICASSO |T§ as well as 
(somewhat model dependent) bounds from Super-Kamiokande [[l!| searches for neutrinos 
from DM annihilations inside the Sun apply [|ll], [17|. Inelastic scattering of a DM particle to 
a nearly degenerate excited state has been proposed in [fT5[| , see ITDL 12^, Q for recent 



analyses, though also in this case tight constraints apply, in particular from CRESST-II [24 



and ZEPLIN-II p5 |. Other proposals include mirror world DM or DM with electric or 



magnetic dipole moments P7| . 

In this work we consider the hypothesis that the DM sector has no direct couplings to 
quarks, only to leptons, in particular the electrons. While electronic events will contribute to 
the scintillation light signal in DAMA, most other DM experiments like CDMS or XENON 
reject pure electron events by aiming at a (close to) background free search for nuclear 
recoils. DM scattering off electrons at rest cannot provide enough energy to be seen in a 
detector, however, exploiting the tail of the momentum distribution of electrons bound in 



an atom may lead to a scintillation light signal in DAMA of order few keV [28]. The signal 



in direct detection experiments from DM-electron scattering has been considered recently 
also in ref . [p9] . An affinity of DM to leptons might also be motivated by recent cosmic ray 
anomalies |30], [31|, [32| observed in electrons/positrons, but not in anti-protons. A simple 
model for "leptophilic" DM has been presented in ref. [33L see in this context for example 
also |54j ^] and references therein. 



In the following we will use effective field theory to perform a model independent analy- 
sis. We will consider all possible Lorentz structures for the effective DM-electron interaction 
and show that in many cases a DM-quark interaction is induced at one or two-loop level by 
photon exchange. In these cases the loop induced DM-nucleon scattering always dominates, 
since the DM-electron scattering cross section is suppressed by the momentum wave func- 
tion. This reintroduces the tension between DAMA and other searches. We identify only 
one possible Lorentz structure, the axial vector type coupling, where DM-electron scattering 
dominates and the scattering cross section is not additionally suppressed by small quanti- 
ties. Taking special care of the kinematics in the DM scattering off bound electrons, we 
show that in this case the fit to the DAMA event spectrum is very bad. Super-Kamiokande 
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FIG. 1: Example for generating an effective local DM-electron interaction vertex (right diagram) 
as used in our analysis by the exchange of a heavy intermediate particle <f> (left diagram). 

constraints on neutrinos from DM annihilations inside the Sun also disfavor this possibility. 
Our results thus suggest that leptonically interacting DM is not a viable explanation of the 
DAMA annual modulation signal. 

The plan of the paper is as follows. In sec. [TT] we introduce the effective Lagrangian for 
DM-lepton interactions, discuss three possible experimental signatures of leptophilic DM, 
and estimate their relative sizes. In sec. [Ill] we discuss all the possible Lorentz structures and 
the implications for tree-level interactions with electrons as well as loop induced DM-nucleus 



interactions. In sec. |TV| the event rates in direct detection experiments are calculated, while 
sec. contains the numerical results for two representative examples, namely vector like 
couplings, where the count rate is dominated by loop induced WIMP-nucleus scattering 
(sec. [V A|) , and axial vector coupling, where no loop contribution is present and WIMP- 
electron scattering off bound electrons dominates (sec. |V Fj| ). In sec. |VT] we show that the cross 
sections required in the axial vector case are ruled out by Super-Kamiokande constraints 
assuming that DM annihilations in the Sun produce neutrinos; we point out the importance 
of the non-zero temperature of the electrons in the Sun. We summarize our results in 
sec. |VII|. Technical details and supplementary information is given in appendices |A|, [B], 0, 
0- 



II. LEPTONICALLY INTERACTING DARK MATTER 



A. Effective dark matter interactions 



The goal of our study is a model independent analysis under the assumption that the DM 
particle x couples directly only to leptons but not to quarks. The appropriate description 
is in terms of effective interactions. Let us first focus on the case of fermionic DM. The 
most general dimension six four-Fermi effective interactions are then, shown pictorially also 
in fig. [I] (right diagram), 

C eS = J2G(xT x x)(im with G=±, (1) 
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where A is the cut-off scale for the effective field theory description, while the sum is over 
different Lorentz structures. A complete set consists of scalar (S), pseudo-scalar (P), vector 
(V), axial- vector (A), tensor (T), and axial-tensor (AT) currents. The four- Fermi operators 
can thus be classified to be of 

scalar-type: T x = c^ + 20^.75, Ti = c e s + icpjs, 

vector-type: T£ = (c v + ^75)7", T tll = (c e v + ^75)7^ (2) 

tensor-type: = (c T + ic AT ^f 5 )a^, T^ v = a^, 

where a^ u = §[7^, ^u]- 1 If DM is a Majorana particle, vector and tensor like interactions are 
forbidden, i.e., Cy = c£ = c\ T = 0. 

In our work we do not rely on any specific realization of the effective interaction. The 
simplest example would just be assuming that the interaction is induced by the exchange 
of an intermediate particle whose mass is much larger than the recoil momenta that are of 
order a few MeV. The intermediate particle can then be integrated out leaving an effective 
point interaction. Let us look at the x-lepton interaction mediated by a scalar field 0, shown 
in fig. |l|. It gives an amplitude 

^s(«x M x) — 2 rr—ids^eut) — ► i (u x u x )(u t u f ) , (3) 

q z — + te 

where on the right-hand side we have neglected the momentum transfer q 2 = (p' —p) 2 <C m|. 
The same amplitude is obtained from a local operator (xx) (^0 with a Wilson coefficient 
dsds/ 171 ^ 11 t ne notation used in eqs. [j], || we have dg = g$, c e s = g e s , A = m^). 

In the case of scalar DM, at lowest order there is only one dimension five operator. The 
effective Lagrangian is given by 

C cQ = G,(x ] x)W s + id Pl5 )i\ with G s = j. (4) 
B. Signals in direct detection experiments 

Let us now discuss the signals that arise when a "leptophilic" DM particle interacts in a 
detector. One can distinguish the following types of signals (see also E9|): 



1. WIMP-electron scattering (WES): The whole recoil is absorbed by the electron that 
is then kicked out of the atom to which it was bound. 

2. WIMP-atom scattering (WAS): The electron on which the DM particle scatters re- 
mains bound and the recoil is taken up by the whole atom. The process can either 
be elastic (el- WAS) in which case the electron wave function remains the same, or 
inelastic (ie-WAS), in which case the electron is excited to an outer shell. 



1 The relation a^j 5 = a a p implies that the AT® AT coupling is equivalent to T®T, and T®AT = 

AT&T. 
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FIG. 2: DM-nucleus interaction induced by a charged lepton loop and photon exchange at 1-loop 
(top) and 2-loop (bottom). 



3. Loop induced WIMP-nucleus scattering (WNS): Although per assumption DM cou- 
ples only to leptons at tree level, an interaction with quarks is induced at loop level, 
by coupling a photon to virtual leptons, see fig. This will lead to scattering of the 
DM particle off nuclei. 

WES produces a prompt electron and possibly additional Auger electrons or X-rays. This 
leads to a signal in scintillation detectors such as DAMA, but is rejected in nuclear recoil 
experiments like CDMS and XENON. If WES was the dominant mechanism, it might be 
possible to explain both the DAMA signal and the absence of the signal in CDMS and 
XENON. In the other two cases, WAS and WNS, the signal consists of a scattered nucleus 
and shows up in all direct detection experiments searching for DM nuclear recoils. If WAS 
or WNS was the dominant signal, then the leptophilic nature of DM would not help to 
resolve the tension between DAMA and the remaining experiments. In the following we 
will first give rough estimates for the relative sizes of these three signal types for different 
DM-lepton effective interactions, while giving a more detailed calculation of the event rates 
later in sec. [TV] and appendix |B|. 

The event rate in direct detection experiments is proportional to the differential cross 
section da/dEa, where 



E v - E' 



(5) 



is the energy deposited by the WIMP in the detector. The DAMA annual modulation signal 
is observed at Ed ~ 3 keV. Also for other direct detection experiments typical values are in 
the few to tens of keV range. As we will see in sec. [TV] and appendix [B|, just from kinematics 



5 




FIG. 3: Momentum space wave functions of iodine and sodium. Solid colored curves correspond 
to shells that contribute to WES in DAMA, while thin light curves are for shells that are not 
accessible in DAMA. 

the cross section is proportional to 

(1(T (xG 2 m e (G 2 m N ) for WES (WAS, WNS) , (6) 



dE d 

where G is denned in eq. |]and m e (mjv) is the electron (nucleus) mass. This suppresses the 
WES induced event rate by a factor m e /mN with respect to WAS and WNS. 

In order for WES to deposit ~ keV energy in the detector, the electron that a WIMP 
scatters off has to have quite a high momentum. Indeed, the maximal detectable energy 
from DM scattering on electrons at rest is 2m e v 2 , with typical DM velocities of v ~ 10~ 3 c. 
Hence, the maximal detectable energy is of order eV, far too low to be relevant for the 
DAMA signal at few keV. Therefore, one has to explore the scattering off bound electrons 
with non- negligible momentum f28fl . In this case the energy transfer to the detector is 
E d ~ 0(pv), and an electron momentum p ~ MeV is required to obtain E d ~ keV. Since 
electrons are bound in the atom, there is a finite yet small probability that it carries such 
high momentum. The detailed calculations below will show that the suppression factor from 
the wave function is given by the expression 

e W ES = ^2m e {E d - E B ) (21 + 1) J ^ \ X m(p)\ 2 ~ HT 8 . (7) 

The integral is over MeV momenta, while Xni{p) is the momentum wave function of the shell 
nl with the binding energy E B . Some wave functions are shown in fig. fj], which we have 
used to obtain the numerical estimate for eyvES ~ 10~ 6 given above. 

Similarly, WAS is also suppressed by the overlap of atomic wave functions of the initial 
and final states of the electron: 

ewAS = \{n'l'm'\e l ^-^\nlm)\ 2 ~ 1CT 19 , (8) 

where the numerical estimate follows from fig. |10| in appendix for a momentum transfer 
of |k - k'| = ^2m N (E d - 5E B ) ~ 10 MeV. 
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Loop induced WNS does not suffer from any wave function suppression, but instead 
carries a loop factor. At 1-loop the suppression is of order (a em Z / V) 2 , with Z being the 
charge number of the nucleus. Combining this with eqs. || [7|, [8|, we obtain the following rough 
estimate for the ratios of WAS, WES and WNS induced event rates (neglecting order-one 
factors but also possible different v dependences): 

R WAS : R WES : R WNS ~ e W AS : e WES — : ( —) ' ~ HT 17 : 10~ 10 : 1 , (9) 

m N V 7T J 

where in the last step we used rriN = 100 GeV and Z = 53. 

We conclude that whenever a loop induced cross section is present it will dominate the 
rate in direct detection experiments. This holds for 1-loop as well as 2-loop cross sections, 
since the latter will be suppressed by another factor (a em Z/7r) 2 ~ 5 x 10~ 6 Z 2 relative to 
1-loop, and hence they are still much larger than the WES contribution. As we will show in 
the next section for axial vector like DM-lepton coupling no loop will be induced. Therefore, 
in this case the signal in DAMA will be dominated by WES. Then, WAS is still irrelevant 
for DAMA, but since WES will not contribute to the rate in CDMS and XENON, WAS 
might in principle lead to a signal in those experiments. 

III. DM ELECTRON SCATTERING VERSUS LOOP INDUCED NUCLEAR RE- 
COIL 

Having estimated the strong hierarchy among the WAS, WES, and WNS signals in the 
previous section, we now discuss which type of signal is present for a given Lorentz structure 
of the effective DM-lepton vertex, eqs. ^ f|. 

A. DM scattering on electrons 

Let us start by investigating DM scattering on electrons, relevant for WES and WAS. To 
simplify the discussion we consider x scattering on electrons at rest. This will enable us to 
see for which types of Lorentz structures in the effective DM-lepton Lagrangian, eq. [I], this 
interaction is relevant. We defer the complications introduced by the fact that electrons are 
actually bound in atoms to sec. [TV] and appendix |B|. 

We consider a DM particle \ of mass m x scattering elastically on a free electron at 
rest, assuming that all the particles are non-relativistic. The scattering cross sections for 
fermionic DM are then: 



scalar- type: a = o\ 
vector-type: a = o\ 




tensor-type: a = {12c 2 -. + 6c 2 AT v 2 } . (12) 

In the above expressions there are two suppression factors, the DM velocity in our halo 
v ~ 10 _3 c and the ratio m e /m x . The cross section for each Lorentz structure is given to 
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TABLE I: Scattering cross section suppression by small parameters for DM-electron scattering and loop 
induced DM-nucleon scattering for all possible Lorentz structures. Here, v ~ 10~ 3 is the DM velocity, 
r e = m e /m x , and qi = mi/rriN (£ = e, fx, r). The reference cross sections er° e , c° e 5 , crLy, c^jv 5 are defined 
in eqs. O, f|l 26. The couplings c x ,c e ,d have been set to one. The entries for \N — > are orders of 



magnitude estimates. 



leading order in these expansion parameters. Up to the velocity or electron mass suppression 
the typical size of the scattering cross section is 



-I - ^ = 4 - 3.1 x ID- 9 cm 2 ( " . (13) 



For scalar DM the x e scattering cross section is induced by the dimension 5 operator, 
eq. |, giving 

° f 4 + fvA, (14) 



2 



with 

a o = Glut _ J_ml _ 7 7 10 - 42 cm 2 / A V 2 / m x \ — 2 

a * e > 5 " 4tt m 2 " 4vrA 2 m 2 " 7,7 * iU Cm ^lOGevj UoOGeW ' 1 &j 

Compared to fermionic DM two powers of A are replaced by m x which typically is larger 
than A. The scalar DM scattering cross section is thus further suppressed compared to the 
fermionic case for given A. The results of eqs. 1C-12 and 14 are summarized in the middle 
column of Tab. |. 
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B. Loop induced DM nucleus interactions 



We have assumed that DM is leptophilic, so that at scale A only operators connecting 
DM to leptons, eqs. [I], @, f|, are generated. However, even under this assumption, at loop 
level one does induce model independently also couplings to quarks from photon exchange 
between virtual leptons and the quarks. The diagrams that can arise at one and two-loop 
order are shown in fig. |2|. 2 The lepton running in the loop can be either an electron or any 
other charged lepton to which the DM couples. 

The one loop contribution involves the integral over loop momenta of the form 



d 4 q 
(4vr) 4 



Tr 



j' + m l h 
q — mi 



m e 



m-„ 



(16) 



with q' = k — k' + q and k, k! the incoming momenta as denoted in fig. ^] and Ti the Dirac 
structures given in eqs. |2j, f|. The one loop contribution is nonzero only for vector and 
tensor lepton currents, 1^ = 7^, a^ v . For the scalar lepton current, 1^ = 1, the loop integral 
vanishes, reflecting the fact that one cannot couple a scalar current to a vector current. The 
DM-quark interaction is then induced at two-loops through the diagrams shown in fig. |2|. 
In contrast for pseudo-scalar and axial vector lepton currents, Ti = 75,7^75, the diagrams 
vanish to all loop orders. One insertion of 75 gives either zero or a fully anti-symmetric 
tensor e a ^ UfJ- . Since there are only three independent momenta in a 2 — > 2 process, two 
indices need to be contracted with the same momentum, yielding zero. 

The calculation of the 1-loop and 2-loop cross sections for scattering of DM on a nucleus 
is relegated to appendix [A} There we give the full 1-loop expressions, whereas here we collect 
the main results in the "leading log" approximation, neglecting the remaining logarithmic 
dependence on momentum transfer. The approximate 2-loop results are obtained in the 
limit of heavy leptons. Expanding also in the x velocity v ~ 10 -3 to first non-zero order, 
the differential cross sections da/dE d are 



vector type: 

tensor type: 
scalar type: 
scalar DM: 



da 



da x N 
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, ( m t 
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21 
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(17) 
(18) 

F{q)\ (19) 
(20) 



where the common 1-loop cross section prefactor is 



dE d 



m N 
2nv 2 



G 



71 



(21) 



2 Similar diagrams with a photon replaced by a Z° or a Higgs boson are power suppressed by (k—k') 2 /M^ „ 
and thus negligible. 
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and da l N 5 / dEd is given by the same expression with G — > G§/{2m x ). Here and Z are 
the nucleus mass and charge, respectively, while //jv = t^n^xI \ m N + w x ) is the reduced 
mass of the two-particle system. The two small parameters are the x velocity v ~ 10 -3 
and the velocity of the recoiled nucleus, Vd = \/2Ed/rriN. The kinetic recoil energy of the 
nucleus Ed in the \N — > xN scattering, cf. eq. |5], has a size Ed ~ keV. 

In the calculations we set \l — A, since this is the scale at which the Wilson coefficient 
G is generated. The form factors F(q) and F(q) account for the nuclear structure. For 
the form factor F(q) entering the one loop induced scattering cross section we use 
K2s2//2 [sin(/€r) 
(and q 2 



36 



3e- 



F{q) 
k = \J1m^E { 



■Kr cos(kt)] / (kt) 3 , with s = 1 fm, r = \J R? — 5s 2 , R = 1.2A 1 / 3 fm, 
^ — k 2 ). The form factor F(q) entering the 2-loop expressions 
accounts for nuclear structure in the case of two-photon exchange. Its precise form is not 
needed in the subsequent analysis, though. The two-loop scalar type differential cross section 
in eq. 19 was calculated integrating out first the leptons assuming they are heavy. This is 
an appropriate limit for muon and tau intermediate states, where m^m T ^> k, while for 
electrons m e ~ k and the expression for the cross section is only approximate, see appendix [S] 
for details. 

For easier comparison with the previous subsection we also quote the results for the total 
xN — ■> xiV cross sections, integrated over the recoil energy Ed- For simplicity we neglect 
the dependence on the nuclear form factors and set F(q) = F(q) = 1 for this comparison, 
giving 



tensor type: a —a^ 



log 




(?)] 




log 


(?)] 





vector type: a =a l N log {^j -{(c^4) 2 + {c x A c l v ) V 

{4 + 4 T -^}log( 
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12/ V m£ 
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(22) 
(23) 
(24) 
(25) 



where a l N is the integral of the differential cross section of eq. El 



7T V 7T / 



x 10 -32 cm 2 



A \ 4 (_Vn_\ 2 ( ^_ 
J VlOGeV/' 1 53 



10 GcV 



(26) 



and cj^g is obtained from the above expression with G — ► G§/ (2m x ). The above results are 
summarized in table |, facilitating comparison with x scattering on free electrons. In table | 
we took yUAr ~ m/vr ~ m x , while the scaling for other values of nucleon and DM masses is 
easy to obtain from above results. 



C. Discussion of Lorentz structures 



In sec. [II B| we have estimated a strong hierarchy between the three types of signals as 
r was < ^wes < ^wns^ gee eq |_ Thege results imply that w henever WNS at 1-loop or 

2-loop is generated it dominates the event rate in direct detection experiments. The Lorentz 
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structures for which this situation applies can be read off from table |. To be specific we 
will use as a representative example of this class the V ® V coupling in the rest of this 
paper. From the table we also see that there is one case — the A ® A coupling — where 
no \N scattering is induced at loop level and moreover the WIMP-electron cross section 
is not additionally v and/or m e /m x suppressed. Hence, we chose the A ® A coupling as 
our second representative example to quantitatively discuss the case of a WES dominated 
event rate. The results from these two examples can be qualitatively extrapolated to other 
Lorentz structures using table | 

As we will see in the following, the \ e ~^ X e cross section in the A ® A case has to be 
very large (corresponding to A ~ (9(100 MeV)) in order to be relevant for DAMA. For the 
cases in table | where <r° e is further suppressed by small numbers, like for example S Cg> P or 
P ®P the scale A would have to be even lower, so that the effective field theory description 
would break down. 

Finally, let us mention the tensor coupling T <g> T, where the 1-loop cross section is 
suppressed by mj/m 2 N , while x e scattering is enhanced by a factor 12. If DM couples only 
to the electron and not to p and r the suppression of the loop is of order m 2 e /m 2 N ~ 10~ 10 , 
and hence, WES and WNS rates can be of comparable size. However, in general one expects 
also a coupling to the p and r leptons. To be specific, in our numerical analysis of V ® V 
and A <g> A cases we will assume equal couplings to all three leptons. For the tensor case the 
same choice would mean that WNS dominates. 



IV. EVENT RATES 

In this section we provide the event rates in direct detection experiments. For WES and 
WAS we assume A ® A coupling and for WNS we take V <E> V. These rates will be used for 
the numerical fits to DAMA, CDMS, and XENON data in the following. As argued above, 
the A <S> A and V <E> V cases are representative enough to cover qualitatively all possible 
Lorentz structures. 

The differential counting rate in a direct DM detection experiment (in units of counts 
per energy per kg detector mass per day) is given by 

dR po rj f 3 da 



d * v — «/ (v), (27) 
dE d m x p det J dE d 

where E d = E x — E' x is the energy deposited in the detector, p Q is the local DM density 
(which we take to be 0.3 GeVcm -3 ), 77 is the number density of target particles, and pdet is 
the mass density of the detector. If the target contains different elements (like in the case 
of the DAMA Nal crystals), the sum over the corresponding counting rates is implied. 

In eq. |27|, /©(v) is the local WIMP velocity distribution in the rest frame of the detector, 
normalized according to J d 3 v /©(v) = 1. It follows from the DM velocity distribution in the 
rest frame of the galaxy, / ga i(v), by a Galilean transformation with the velocity of the Sun in 
the galaxy and the motion of the Earth around the Sun. For / ga i(v) = f g3 \(v) we assume the 
conventional Maxwellian distribution with v = 220 km s" 1 and a cut-off due to the escape 
velocity from the galaxy of f 0S c = 650 kms -1 : f gSL \(v) oc exp(— v 2 /v 2 ) — exp(t>^ sc /'u 2 ) for 
v < v esc and zero for v > v csc . We have checked that the precise value of the escape velocity 
has a negligible impact on our results. 
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A scattered nucleus does not deposit all its energy in the form of scintillation light. This 
effect is taken into account by the so-called quenching factors, which are gNa — 0.3 for sodium 
recoils and q\ ~ 0.085 for iodine recoils |37| in DAMA. In refs. [39[ it has been pointed out 
that the so-called channeling effect could be relevant, implying that for a certain fraction of 
events no quenching would occur due to the special orientation of the recoil with respect to 
the crystal. So far this effect has not been confirmed experimentally in the relevant energy 
range. Indeed, the results of ref. [EDI do not indicate the presence of any variation of the 



count rate for special crystal directions. Quenching and channeling is relevant in DAMA 
in the cases of WAS and WNS , while the scattered electrons in the case of WES produce 
unquenched scintillation light. In our fit to DAMA data for WNS we do include channeling 
following ref. |39| (similar as in [0), but we comment also on the case when no channeling 
occurs. 



A. WIMP electron scattering 

To obtain an expression for the event rate in the case of WES it is necessary to take 
into account the fact that electrons are bound to the atoms. The kinematics of scattering 
off bound electrons has some important differences compared to scattering off free particles. 
The bound electron does not obey the free-particle dispersion relation E^s = p 2 + m 2 . 
Instead it has a fixed energy E e = m e — Eb, determined by the binding energy of the atomic 
shell, Eb > 0, whereas its momentum p follows a distribution which is given by the square of 
the Fourier transform of the bound state wave function corresponding to that shell. Energy 
conservation reads in this case E x + m e — Eb = E' + E' e , or 

E' e = m e + E d -E B . (28) 

After some tedious algebra one arrives at the following expression for Ed- 

2 

E d ^ — pvcosO, (29) 

2m x 

where 3 cos 9 = kp/kp and we used the approximation E d <C m e < E e <C m x and v ~ 10~ 3 . 
We see that to obtain detectable energies relevant for DAMA (Ed of few keV), electron 
momenta of order MeV are required. 

In appendix |B| we give the details on the calculation of the scattering cross section and 
count rate in the case of WES, taking into account the peculiarities of scattering on bound 
electrons. Here we only report the final results. Assuming the axial vector Dirac structure, 
T x = T e = A, as motivated above, the count rate is (we also set = c e A = 1 for simplicity): 



dR WES 3p m e G 2 ^ r — — " f dpp 



dE d 4vr(m I + m Na )ra x , 



J2^rn e (E d -E B , nl )(2l + l) J ' ^ \ Xn i( P )\ 2 1(v™) 

(30) 



We always denote the DM momentum with k and the electron (or nucleus) momentum with p, see fig. 0. 
Bold symbols refer to 3- vectors and k = |k|, and similar for p. 
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Here Xniip) is the momentum wave function of the electron, and the function I(v m i n ) is 

I(v min ) = [ d 3 v^- 9(v - v min ) , (31) 



while the minimal velocity required to give detectable energy Ed follows from eq. 



WES + _JL_ I ':?') ) 

^ min p + 2m x • {6I) 



For m x > 10 GeV and p of order MeV the first term dominates. 



The sum in eq. ^0] is over the atomic shells of both iodine and sodium with quantum 
numbers nl, and Eb^i is the corresponding binding energy. The electron can only be kicked 
out of its atomic shell if its binding energy is smaller than the total energy deposited in the 
detector (cf. eq. |28| ): 

Ed > E B ^ n i . (33) 

Only the shells satisfying this requirement can contribute to the event rate in eq. The 
momentum wave function Xni{p) is defined in eq. B3] in appendix O. Technical details on 



how we implement the wave function numerically are given in appendix [C[ the results for the 
iodine and sodium wave functions are shown in fig. |3|. We see that the dominant contribution 
to WES scattering in DAMA comes from the inner s-shells of iodine because these are largest 
at high p. Electrons from the Is, 2s, 2p shells — depicted as thin light curves in fig. |3] - 
do not contribute to the DAMA signal region of Ed ~ 2 — 4 keV since the binding energies 
are too large, respectively 33.2 keV, 5.2 keV, and 4.7 keV [fl~I}| . The shell dominating the 
signal in the 2-4 keV region is the 3s shell of iodine, with a binding energy of about 1 keV. 
Apparently this has been overlooked in ref . [ 2£| , while it has important consequences on the 



size of the needed cross section, see discussion in sec. IV B 



B. WIMP atom scattering 

Let us consider now the case when the electron on which the DM particle scatters remains 
bound and the recoil is taken up by the whole atom. According to the coordinate space 
Feynman rules, the matrix element for WAS for an electron in atomic shell nlm in the initial 
state and n'l'w! in the final state is given by 

M r ^ M = G J A4, m ,(x)i im (x)e- ip ' x e ! ( k - k ') x <rX«;'r e X. (34) 

Here, k and k' are the initial and final momenta of the WIMP, and p' is the average 
momentum of the electron in the final state resulting from the motion of the whole atom 
after the scattering. Since most of the recoil momentum is carried by the nucleus, |p'| is 
smaller than |k — k'| by a factor of m e /mjy, and can therefore be neglected. The coordinate 
space wave function of the electron in the state with orbital quantum numbers nlm is denoted 
by i>nim(' x ) ■ Again we specialize to the case of axial vector coupling, = T^f = and 
set c\ = c e A = 1. We use non-relativistic spinors, which is certainly justified for u r x and u r x \ 
and also for u s e except, perhaps, for electrons from the Is shell of iodine. In this last case, 
relativistic corrections are of order 20%. 
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Let us first consider the case when the electron remains in its state, and hence the scat- 
tering on the atom is elastic (el- WAS). Then we have s — s' and nlm = n'l'm' . Furthermore, 
we have to sum coherently over all shells and electron spins, since it is impossible in principle 
to identify on which electron the WIMP has scattered. It turns out that for the axial vector 
case the spin sum vf^^ul vanishes. This can be verified by using explicit expressions 
for the spinors u s e , and follows from the fact that the different sign due to 75 of right-handed 
and left-handed components of the electron cancel each other in case of a coherent sum 
over spins. 4 The elastic scattering may be relevant for other Lorentz structures where this 
cancellation does not occur. However, in sec. [Ill C| we have argued that the only case of 



practical relevance is the axial coupling, and therefore we will not consider el- WAS further. 

We are left now with the case where the electron is excited to an outer free shell which 
corresponds to inelastic WIMP-atom scattering (ie-WAS). In this case the sum over all 
occupied electron states nlm, over all unoccupied states n'l'm', and over WIMP and electron 
spins has to be incoherent because one can distinguish in principle different initial and final 
states, e.g. by x-ray spectroscopy. The differential cross section in this case is obtained as 



da WAS m N \M\ 2 

(35) 



dEd 32Tcm 2 m'iv 2 



Plugging in the matrix element from eq. 133 we get 

da ie ~ WAS 3m N G 2 



with 



2?n , £ £ \(n'l' m y^\nlm)\* (36) 

nlm n'l'm' 



(n'l'm'\e^-^\nlm) = J d 3 x<, IW (x)Vw(x) e^'^ . (37) 



The expression for the counting rate is obtained from eq. [27], 

A Die- WAS o „ ri2 

dK^_ = m N S P0 G y y | (nYm y(k-k')x| n , m) |2 /( ^WA S) ( g g) 

dE d m l + m N ,27im x ^^^ 1 K } 

^ nlm n'l'm' 



with N = I, Na. The function / is defined in eq. |31], and the minimal velocity required 
to give detectable energy Ed follows from the kinematics implied by energy conservation, 
Ed = E x — E' = 5E B + m N v 2 N /2, and momentum conservation, k = k' + m N v N : 

^ie-WAS = E d{m x + m N ) - m N 5E B 
m x ^2m N (E d - 5E B ) 

where 6E B is the difference of the binding energies of the initial and final shells: 6E B = 
E Bt nim ~ E Bt n'i'm'- Details on how we calculate the matrix elements involving the wave 
function in eq. |38] are given in appendix 0. 



4 This argument will not hold if an unpaired valence electron is available so that we cannot sum over spins. 
However, most chemically bound systems are formed in such a way that this does not happen. Even 
in this case el- WAS would be suppressed since scattering on outer electrons is highly suppressed by the 
smallness of the binding energy of these electrons compared to the transferred momentum. 
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C. Loop induced WIMP nucleus scattering 



The event rate for loop induced DM-nucleus scattering is easy to obtain from the differ- 
ential cross sections daj^ jdE d in eqs. |17H2~0| and the general expression for the counting rate 
eq. 

' d<j jv 



jpWNS 

dE d 



Po 



1 



,WNS\ 



m x mi + m Na \dE d ~ )' v ~ mm ' ' ^ 
while the minimal velocity to produce a detectable 



with 



The function I is defined in eq. pi , 

energy Ed is given for WIMP-nucleus elastic scattering by v^^ s 
jj N = m x m N I (m x + m N ). 

We now specialize to the V <S> V case. The event rate depends on the x mass and the 



^E d m N /2/j 2 N 



c 



V 



coupling constant of the effective operator G (we set 
comparison with previous works, it is useful to trade G for the total x e 

rO 



1 from now on). For easier 



Xe cross section 



(7° = G 2 m 2 ./7i, eq. |13[ For the V <8> V case we thus have 



da 



dE d 



N 2 
V 



^(-)W 



log -f 



(41) 



to be inserted in eq. As discussed in appendix [S] this leading log approximation is quite 
accurate. Nevertheless, in the numerical calculations we use the full expressions given in 
appendix 0. We set \i = 10 GeV, since this corresponds roughly to the scale A, where our 
effective theory is defined. Furthermore, we assume (somewhat arbitrarily) equal couplings 
to all three leptons. The logarithm in eq. f|l] implies then a relative contribution of e : fi : 
r ~ 30 : 7 : 1. Note that the rate is dominated by the contribution from the electron in 
the loop assuming equal couplings at the scale A ~ 10 GeV. Therefore, our results are 
conservative, in the sense that per assumption DM has to couple to the electron. 



V. FIT RESULTS FOR DAMA, CDMS, AND XENON10 

A. Vector like interactions and loop induced WIMP nucleon scattering 



The event rate in the case of the loop induced WNS, eq. |40], is very similar to the 
corresponding expression for usual spin-independent elastic WIMP scattering on the nucleus 
(see, e.g., refs. [|TI| (TTJ, P5[), denoted as "standard case" in the following. The main difference 
is the replacement of the atomic mass number A by the charge number Z (and an additional 
logarithmic E d dependence beyond the leading log approximation). Therefore we expect 
that the fit of DAMA and the compatibility to CDMS and XENON will be very similar to 



the standard case, see e.g., refs. M, 10, \L \ 



Our numerical analysis of DAMA, CDMS, and XENON data follows closely ref. JTC 
where technical details on the fit can be found. For the DAMA fit we use the spectral data 
on the annual modulation amplitude S m from the threshold of 2 keV up to 8 keV. The data 
above this energy are consistent with no modulation and since our model does not predict 
any features in that region they do not provide an additional constraint on the fit, apart 
from diluting the overall goodness of fit. In addition to the data on the amplitude of the 
modulated count rate we use also the unmodulated event rate as a constraint in the fit. While 
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FIG. 4: Predicted spectrum for the modulated (top) and unmodulated (bottom) event rate in DAMA 
at the best fit point assuming loop induced WIMP-nucleus scattering resulting from vector-like DM-lepton 
couplings. Results are shown for using all data points from 2-8 keV (solid) and for omitting the 1st bin 
(dashed). The parameter values and the x 2_v alues are given in the legend. 



the bulk of these events will come from various unidentified backgrounds, every model has 
to fulfill the constraint of not predicting more unmodulated events than actually observed in 
DAMA. Fig. |] shows the predicted spectrum at the best fit point compared to the DAMA 
data. Note that the error bars on the unmodulated rate (lower panel) are hardly visible for 
most of the data points, as a result of the huge number of events in DAMA. For the analysis 



using data from 2-8 keV the best fit point is at 



12.4 GeV, o\ e 



4.5 x 10" 44 cm 2 , and 



we obtain the excellent fit of x 2 /dof = 9.1/10. If we drop the first data point the fit even 
improves to x 2 /dof = 2.8/9. 

The allowed regions in the plane of DM mass and scattering cross section are shown 
in fig. |5|. For easier comparison with the case of scattering off electrons we parameterize 
the cross section on the vertical axis in terms of erJL , see eq. 41. The spectral data on the 
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FIG. 5: DAMA allowed region at 90% and 3a CL in the case of 1-loop induced WIMP-nucleus scattering 
(V ® V coupling) in the plane of the WIMP mass and the WIMP-electron cross section tr° e = G 2 m1/Tr. 
Regions are shown with and without taking into account the channeling effect. Furthermore, we show the 
bounds at 90% CL from CDMS-II and XENON10. The dashed curves show the 90% CL constraint from 
the Super-Kamiokande limit on neutrinos from the Sun, by assuming annihilation into rf or vv 1 see sec 
for details. 



VI 



modulated signal results in an allowed region for rather small DM masses around m x ~ 
12 GeV. If channeling is not taken into account an allowed region appears at similar DM 
masses but at higher cross sections (due to scattering off sodium J42|). 5 This DAMA allowed 
region has to be compared to the constraints from CDMS and XENON. The compatibility 
with these bounds is similar to the standard case: while marginal compatibility might remain 
there is clearly severe tension between the DAMA signal and the CDMS and XENON 
bounds in this framework. Here we will not elaborate on this question further and refer to 



refs. 0, 0, HI for detailed discussions of the DAMA versus CDMS/XENON compatibility 



in the standard case. 

The main motivation for considering electron interacting DM — namely avoiding the con- 
straints from nuclear recoil experiments, is thus invalidated by the loop induced nucleon 
scattering. We now turn to the axial-axial coupling, where the loop induced scattering can 
be avoided. 



5 In both cases (with and without channeling) there is also a local minimum around a WIMP mass of about 
80 GeV from unchanneled scatterings off iodine. In fig. ^ we show confidence regions defined with respect 
to the global minimum, and this second region does not appear at 3er if channeling is included. 
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FIG. 6: Predicted spectrum for the modulated (top) and unmodulated (bottom) event rate in DAMA at 
the best fit point assuming WIMP-electron scattering resulting from axial vector-like DM-electron couplings. 
Results are shown for using all data points from 2-8 keV (solid) and for omitting the 1st bin (dashed). The 
parameter values and the y 2 -values are given in the legend. 

B. Axial vector like interactions and WIMP electron scattering 



We perform a similar fit to DAMA data as before but using now eq. ^ f° r the event 
rate for WES. The predicted modulated and unmodulated DAMA event rates at the best 
fit in this case are shown in fig. || Using the data from 2-8 keV we obtain a rather bad fit 
to the modulated spectrum with x 2 /dof = 55.9/10, which corresponds to a probability of 
2 x 10~ 8 . The prediction drops too fast with energy in order to provide a satisfactory fit to 
the data. If we omit the first energy bin the fit improves considerably to x 2 /dof = 20.6/9 
corresponding to a probability of 1.4%. We find, however, that the parameter values from 
this fit predict a very sharp rise for the spectrum of the unmodulated event rate in DAMA, 
see lower panel of fig. || In the fit we have required that the unmodulated prediction stays 
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FIG. 7: DAMA allowed region at 90% and 3er CL in the case of WIMP-electron scattering (axial-axial 
coupling) in the plane of the WIMP mass and the WIMP-electron cross section <j° e = G 2 m} e j-K. We 
stress that these regions are obtained with respect to the best fit point, which by itself does not provide a 
satisfactory fit to DAMA modulated and unmodulated spectral data, see fig. [| We show also the bounds at 
90% CL from CDMS-II and XENON10 from inelastic WIMP-atom scattering. The dashed curves show the 
90% CL constraint from the Super-Kamiokande limit on neutrinos from the Sun, by assuming annihilation 
into rf or vv, see sec. |Vl] for details. 

below the observed rate within the analysis window down to 2 keV. However, DAMA shows 
also some data points for the unmodulated rate below 2 keV, which are not compatible with 
the predicted rate. While it is not possible to use data below 2 keV for the modulation, it 
seems likely that they rule out models predicting more events than observed. The WES fit 
shown as dashed curve in fig. |6| predicts more than a factor 3 more events than observed 
in the first two bins below 2 keV, where error bars are still very small. We conclude that 
WES has severe problems to explain the spectral shapes of the modulated and unmodulated 
components of DAMA data. 

If we ignore the problems of the spectral fit and despite the low goodness-of-fit consider 
"allowed regions" in the plane of WIMP mass and cross section relative to the best fit point 
we obtain the results shown in fig. [7[ We observe that very large cross sections are required: 

*~ 1 ^ <m ' x (i55&v)- < 42 > 

where the linear dependence on m x holds for m x > 10 GeV. The vastly different best fit 
cross sections for WNS and WES follow from the discussion in sec. [II B| where we estimated 
the relative size of the corresponding counting rates, see eq. 0. Here we do not explore other 



phenomenological consequences of such a large cross section. Just note from eq. [13] that we 



can realize a cross section of this order of magnitude only with a relatively low scale for the 
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new physics of A < 0.1 GeV, where we have assumed for the couplings ~ c A ~ 1. If 
in a particular model constraints on the coupling constant c e A apply (see, for example [[J3|), 
A has to be accordingly smaller. Note that the momentum transfer for WES is given by 
the momentum of the bound electron, which has to be of order MeV. This provides a lower 
bound on A in order to describe the interaction by using the effective theory. 

The cross section from eq. ^2] is about 5 orders of magnitude larger than the result 
obtained in ref. P5[ , which finds pb-size cr° e at m x ~ 100 GeV. Let us comment, therefore, 
on the differences of our analysis to the one from |28| . Apparently the main difference is that 



we take into account the special kinematics related to the scattering off bound electrons, 



whereas in ref. [28] electrons are treated as effectively free with a momentum distribution 



obtained from the wave function. Our calculation outlined in sec. |IV A| and appendix [B| 
leads to several suppression factors in the WES event rate with respect to the expression 
used in [28]. Our minimal velocity v^^ s from eq. p2| is roughly a factor two larger than the 



one used in p8| requiring to go further out in the tail of the WIMP velocity distribution. 
Furthermore, the condition eq. [33|, > Eg, which prevents the contribution of the inner 
shells of iodine, has not been imposed in |28||. From fig. 13 we see that at p ~ 1 MeV the 



wave function of the iodine Is shell is about 2 orders of magnitude larger than the one of the 
3s shell, which actually gives the first relevant contribution after requiring that the binding 
energy has to be lower than E^ ~ few keV. 

As mentioned above, in the case under consideration inelastic WAS may contribute to ex- 
periments searching for nuclear recoils. To calculate this effect we would have to perform the 
sum in eq. ESI over all occupied shells nlm and all free shells n'l'm' . It turns out numerically 



that transitions from s-shells to s-shells give the largest contributions, see also appendix y. 
In order to estimate the order of magnitude we have taken into account transitions from the 
Is, 2s, 3s orbitals to the first two free s-orbitals of the germanium and xenon nuclei relevant 
in CDMS and XENON, respectively. In fig. [7| we show the constraints resulting from this 
estimate of the ie-WAS event rate. Numerically these constraints turn out to be very weak 
and the limits are several orders of magnitude above the region indicated by DAMA; the 
good sensitivities of CDMS and XENON to nuclear recoils cannot compensate the large 
suppression of the ie-WAS count rate compared to WES, as estimated in eq. [| 

Although the poor quality of the fit in the case of WES already disfavors this mechanisms 
as an explanation for the DAMA modulation signal, we will show in the next section that 
constraints on neutrinos from DM annihilations inside the Sun are even more stringent 
and exclude the cross sections required for DAMA by many orders of magnitude if DM 
annihilations provide neutrinos in the final states. 



VI. NEUTRINOS FROM DM ANNIHILATIONS INSIDE THE SUN 



Any DM candidate has to fulfill the constraints on the upward through-going muons 
coming from water Cerenkov detectors, like Super-Kamiokande [[H|, and from neutrino 
telescopes |45], f46|. Some recent papers [[17], [47j have discussed the constraints on the 
DAMA region in the framework of standard WIMP-nucleus scattering. Here, we reanalyze 
the bound coming from the Super-Kamiokande experiment in the framework of leptonically 
interacting DM. 
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FIG. 8: WIMP capture rate in the Sun as a function of the WIMP mass assuming DM scattering off 
electrons, hydrogen, and all other nuclei in the Sun, with a scattering cross section of IO -36 cm 2 . The solid 
curves correspond to scattering off particles at zero temperature, whereas the dotted curves show the effect 
of the actual temperature distribution inside the Sun [ fl8[ for electrons and hydrogen. 



One important ingredient for the prediction of the neutrino flux coming from DM an- 
nihilations inside a celestial body is the capture rate C , which is proportional to the DM 
scattering cross section, see, e.g. j56fl . In the calculation of this quantity, usually, the WIMPs 
are assumed to interact with material at zero temperature, neglecting the solar temperature 
of about 1.5 x 10 7 K in the center and 8.1 x 10 4 K at the surface. Although this is a reason- 
able assumption for WIMP candidates interacting with hydrogen and the other nuclei inside 
the Sun, it fails for the case of DM scattering on the free electrons in the Sun. The effect 
of non-zero temperature on the capture rate depends on the ratio of the thermal velocity 
of the target to the WIMP velocity. The thermal kinetic energy fceT is independent of the 
mass, but the thermal velocity is larger by a factor ^Jm p /m e ~ 45 for electrons compared 
to hydrogen. 

We calculate the rate for WIMP capture by a body at finite temperature following ref. p5 



where the expression given there has to be extended to include the motion of the Sun with 
respect to the DM halo. The temperature distribution for the electrons inside the Sun is 
taken from the solar model BS2005-AGS,OP f48fl . Fig. § shows the effect of the non-zero 
temperature on the capture rate for electrons, hydrogen and all other nuclei in the Sun. We 
find that the capture rate on electrons is enhanced by about one order of magnitude, while 
the effect is hardly visible at the scale of the plot for hydrogen. The temperature effect can 
be neglected for scattering off heavier nuclei, which dominates the capture in the case of loop 
induced WNS. In this case one has to include also a suppression due to the nuclear form 
factor. Furthermore, we neglect the effect of WIMP evaporation, important only for DM 



masses lower than 10 GeV [50] and the gravitational effects from planets like Jupiter [51 
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The annihilation rate T is related to the capture rate C© by |36 



C© , i 2 / > ' 



^ = -t^\j)' T =vu^- (43) 

where t® ~ 4.5 Gyr is the age of the Sun and the parameter Ca depends on the WIMP 
annihilation cross section and on the effective volume of the confining region in which the 
DM particles are trapped: C A = (a^v) V 2 /V? with Vj = 6.6 x 10 28 (j m x /10 GeV)~ 3 / 2 cm 3 . 
We denote by (a ann v) the thermally averaged total annihilation cross section times the 
relative velocity, at the present time. Capture and annihilation are in equilibrium if r <C t Q . 
Then the annihilation rate is just half the capture rate and becomes independent of the 
annihilation cross section {a ann v). For our calculations we assume this limit, we comment 
on its validity in appendix [D]. 

The neutrino flux at the detector from the annihilation channel / with branching ratio 
BRj is given by 

dti _ r dNf 

dE v - BR ^lE v > (44) 
with d being the distance between the Earth and the Sun. Here we are interested in annihi- 
lations into leptons. We consider the following four channels: rf , u e u e , u^P^, and v T v T . Note 
that annihilations into electrons do not provide neutrinos, and muons are always stopped 
before decay, giving rise to neutrinos in the MeV energy range which is below the Super- 



Kamiokande threshold [ 5"2"| . In the case of direct neutrino channels, the initial neutrino 
spectrum is simply a Dirac 5 function centered at E v = m x , and we assume a flavor-blind 
branching ratio, i.e., BR UePe = BR^p^ = BR UtPt = 1/3. The results do not depend strongly 
on this assumption, since flavors are mixed due to oscillations. 6 For the rf channel, we use 



the initial neutrino spectrum given in | 53fl 



The neutrino spectrum dN£ /dE u at the detector is calculated considering the effect of 
neutrino oscillation, coherent MSW matter effect, absorption, and regeneration (see e.g., 
|53| , [54}] ) by solving numerically the evolution equations of the neutrino density matrix 
within the Sun. The neutrino oscillation parameters and mass squared differences are fixed 
to the best fit values reported in |55]. We set #13 to zero, avoiding in this way possible Earth 



matter effects. 

The total muon flux is given by the formula, see e.g. ]56 



$ M = / dEj dE v ^ N A R,(E„ E i 

JEt h JE U atj u 



da?, ,„ „ s . , da!, 



d( t>l » r 7, ™th\ 

(45) 

where M p and M n are the fractional number of protons and neutrons at the point of muon 
production. We consider M p — Af n ~ 0.5, since for the through-going muons in the Super- 
Kamiokande detector the interaction can be assumed to occur in standard rock for which 



6 There is some difference of the z/ T j/ r -channel due to r regeneration effects, which are important for high 
energies. Assuming annihilations with branching ratios equal to one for each of the three flavors we 
find that the muon neutrino flux at the Earth is practically the same for all three initial flavors up to 
m x ~ 100 GeV. For m x = 1 TeV the ratio of the muon neutrino fluxes at Earth is roughly 1 : 3.5 : 6.4 for 
annihilations into v e v e w^v^: v T v T . 
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the number of protons is almost equal to the number of neutrons (Z = 11, A = 22). Na 
is Avogadro's number, and the effective number of nucleons per unit volume is given by 
NAProck, with p roc k the density of the material. The muon range is defined as the 
distance traveled by a muon with initial energy and final energy equal to the detector 
energy threshold E*' . 

=M£I|) (46) 

with a ~ 2.2xlO~ 3 GeV/ (gem -2 ) and (3 ~4.4xl0~ 6 / (gem -2 ). The muon energy threshold 
has been fixed to Ejj 1 = 1.6 GeV, corresponding to the 7 m path- length cut applied on upward 
through-going muons in the Super-Kamiokande detector. For the differential cross sections 



daP' n / dEp, we use the analytic expressions for deep inelastic scattering given e.g. in ref. [56 



Super-Kamiokande gives 90% CL limits on the muon flux induced by neutrinos from DM 
annihilations in the Sun for cones of different opening angles around the direction of the 
Sun [16]. To be conservative we use the limit for a cone with half- angle of 20°, which should 
include 90% of all muons at m x ~ 18 GeV [IB|, and a fraction approaching 100% for larger 
WIMP masses. The corresponding limit is $ M < 1.1 x 10 -14 cm" 2 s -1 fl6f . Using this upper 
bound on $^ we obtain via eq. f|5] an upper bound on the DM scattering cross section as a 
function of m x . 

These bounds are shown in figs. |^ and [7] for the case of loop induced WIMP-nucleus 
scattering (WNS) and WIMP-electron scattering (WES), respectively. We show the limit 
for annihilations into rf and vv (assuming equal branchings into the 3 flavors). In the 
case of WNS, annihilations into neutrinos exclude the region compatible with DAMA, while 
annihilations into tau leptons might be marginally consistent with it at 3a. In contrast, in 
the case of WES the neutrino bound excludes the region indicated by DAMA by more than 6 
orders of magnitude. This implies that if DM couples to electrons with a cross section as large 
as indicated by the WES DAMA fit, c.f. eq. W\ DM annihilation into neutrinos must be very 
strongly suppressed. This will be hard to achieve because annihilation into charged leptons 
generates almost model independently also annihilation into neutrinos from W exchange at 
1-loop. Thus annihilation into neutrinos is typically suppressed by a loop factor of O(10~ 4 ) 
compared to annihilation into charged leptons. This rules out all leptophilic DM models 
with dominant direct annihilation into leptons as an explanation of DAMA/LIBRA. 



VII. CONCLUSIONS 

In this study we have considered the hypothesis that DM has tree level couplings only 



to leptons but not to quarks. Such a model has been proposed in ref. [£§] to reconcile the 
DAMA annual modulation signal with constraints from searches for nuclear recoils from DM 
scattering. Our results imply, however, that this is not possible for the following reasons: 

1. By closing the lepton legs to a loop, we obtain a coupling to the charge of the nu- 
cleus by photon exchange. Whenever the Dirac structure of the DM-lepton coupling 
allows such a diagram at 1 or 2-loop WIMP-nucleus scattering will dominate over 
the scattering rate from the direct coupling to electrons, because the latter is highly 
suppressed by the high momentum tail of the bound state wave function. This leads to 
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a situation very similar to the standard WIMP case, implying the well-known tension 
between DAMA and the constraints from CDMS and XENON10, see fig. §. 

2. If the DM-lepton coupling is axial vector like, no loop will be induced and hence the 
scattering proceeds only by the interaction with electrons bound to the atoms of the 
detector. We have performed a careful analysis of this case, taking into account the 
peculiarities of scattering off electrons in bound states. We find that this model is 
strongly disfavored as an explanation of the DAMA signal because 

(a) the predicted spectral shape of the modulated and/or unmodulated signal in 
DAMA provides a very bad fit to the data as shown in fig. |6], and 

(b) the cross section required to explain the DAMA signal is ruled out by Super- 
Kamiokande constraints on neutrinos from DM annihilations in the Sun, see 
fig- 0- 

The arguments 1 and 2(a) are rather model independent, relying only on the presence of 
the effective DM-lepton vertex, while the argument in 2(b) depends on the assumption that 
neutrinos are produced by DM annihilations. Due to SU(2)l gauge symmetry, generically 
one expects that DM will couple to both, charged leptons and neutrinos, which would open 
the annihilation channel into vv. If for some reason DM couples only to charged leptons, DM 
would generically also annihilate into rf, leading again to the neutrino signal. In order to 
evade the Super-Kamiokande constraint one has to forbid the coupling of DM to neutrinos 
and to the tau lepton. Let us mention that the most generic way to avoid coupling to 
neutrinos is the chiral coupling only to right-handed leptons. Note, however, that such a 
chiral V + A coupling involves a vector-like coupling which will induce DM-quark scattering 
via the loop diagram and argument 1 applies. Another way to evade the bound from 
annihilations would be to assume that DM is not self-conjugate and postulate the presence 



of a large x~ X asymmetry in our halo, see e.g., refs. |57], |58|, |59 . 

In conclusion, we have shown that the hypothesis of DM-interactions only with leptons 
does not provide a satisfactory solution to reconcile the DAMA annual modulation signal 
with constraints from other direct detection experiments. 



Note added: After the completion of this work we became aware of Ref. | 60fl , where 
CDMS publishes constraints on electron-like events above 2 keV in their detector. These 
results apply to the case of axial coupling, where scattering off electron dominates. From 
Fig. 6 we find that our fit predicts an unmodulated rate of 0.4 (1.5) events/d/kg/keV at 
2 keV if the lowest energy bin of the modulated rate is (is not) taken into account. CDMS 
observes 1.93 ± 0.24 events/d/kg/keV at 2 keV. Assuming a flat background in the energy 
range of interest an upper limit on a possible signal from DM of 0.5 events/d/kg/keV at 
90% CL is obtained at 2 keV (see Fig. 3 of f60|). The signal in Ge is expected to be 
roughly one order of magnitude smaller than in iodine due to wave function suppression. 



This estimate suggests that the results of Hi do not rule out the DAMA region shown in 



Fig. 7. A more detailed analysis may still be of interest, though. 
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APPENDIX A: LOOP INDUCED DM QUARK INTERACTIONS 

In this appendix we calculate the cross sections for DM-nucleon scattering through the 
loop-induced interactions shown in fig. |2|. The main results were already collected in sub- 
section [111B| in the leading log approximation, while here we give full 1-loop results and 



describe how the approximate 2-loop results were obtained. For calculations we use the 



FeynCalc package [61]. The cross section for scattering of a non-relativistic DM particle \ 



with mass m x on a nucleus at rest carrying charge Z and having a mass is 



\MV- 



dE d 32ixm N m 2 v 2 ' 

A 

with v ~ 10 -3 the x velocity, the kinetic recoil energy of the nucleus and M. the matrix 
element for \N — > xN scattering. 

We start with the vector type interaction between leptons and DM, £g = 
G(xr^x)(£cy7 /i £) with = (c v + 0^75)7^, see eq. |l[ The matrix element for \N — > \N 
scattering generated through the one loop diagram of fig. |2| is then 

i (A2) 
= C^\^){u' x Y^u x )ZF{ q )(u' Nlil u N ) . 

The sum is over the light quarks qi with charges Qi, F(q) is the nuclear form factor defined 
in sec. [llllj| , and Cy (fi) is the 1-loop factor calculated in the MS scheme 



= 2«em Gc ^ 1 dxx ( l _ x)log -V- -_T± ; -I ^ (A3) 



TT 







—x(l — x)q 2 + mj — iO 



where q 2 ~ —k 2 = —Im^Ed is the momentum transfer (in the calculation of the Super- 
Kamiokande bounds we used q 2 ~ — O^m^v 2 )). In the calculation we set \i = A, with 
A ~ 10 GeV, because this is the scale at which the Wilson coefficient G is generated. 7 For 
Tfig 3> k one can neglect the momentum transfer in the above integral, giving an approximate 
expression 

C v L (») = ^G4log(m}/^), (A4) 



7 This choice of [i does not minimize the size of the logarithm in Cy \[i). However, since the expansion 
parameter a om is small this does not invalidate the use of perturbation theory. For a choice of ^< A one 
would need to take into account renormalization group flow and mixing of operators. 
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FIG. 9: The real (solid line) and imaginary (red dashed line) parts of the loop factor ratios rc = /C LL 
as a function of q 2 /m 2 for vector (thick lines) and tensor (thin lines) lepton currents (these two lines overlap 
within the precision that can still be seen on the plot). The left plot is for q 2 negative (space-like momentum 
exchange), the right for q 2 positive (time-like momentum exchange). We take mi = m e and ^ — 10 GeV. 



which is very precise for muon and tau running in the loop. It is quite precise also for 
the electron, even though m e ~ k. The reason is that there is still a hierarchy m e <C 
/i ~ A. Neglecting the difference between m e and k then corresponds to a leading log 
(LL) approximation in the renormalization group running, while the induced error is only 
logarithmic in 1 + 0{n/m e ). 

The \N —> xN differential cross section da/dE^ in the leading log approximation is 
given m eq. P. Multiplying it by \C ( y ] /C V L \ 2 one obtains the full 1-loop prediction. In fig. 
we show the value of Cy /C V L for rri£ = m e and fi = 10 GeV. Note that above the pair 
production threshold, q 2 > Ami, ^ develops an imaginary part, since electrons in the loop 
can go on-shell. The important thing for our purposes, though, is that Cy /C V L is a slowly 
varying function of q 2 and is of (9(1) in the range of space-like q 2 ~ —m 2 of interest to us, so 
that the LL approximation is quite precise. Even so, in the numerical analysis in section 
we use the full 1-loop results. 

Let us next move to the tensor DM-lepton interaction, Ce = G^x^^X)^^) with 
T 1 ^ = (ct + icATlb )u^ u . The matrix element for xiV — > xiV scattering is 

M = C T {p) [u' x Y^u x ) q - 2 (N(p') | Qi (ft7/.ft) \N(p)) 

(A5) 

= Cr(p.) K^X) ^ZF(q) (u' Nltx u N ) , 

with the 1-loop factor in dimensional regularization (as before the pole is to be subtracted 
using MS scheme) 

c (l) = _2aem m ^£ o(g 2 jTO 2 )TO 2) (Ag) 
7T 

Here B (k 2 ,m 2 ,M 2 ) is a two-point scalar Veltman-Passarino function. Explicit expressions 



for it can be found e.g. in [62|. In the leading log approximation the above expression 
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becomes (in the MS scheme) 



C LL (/i ) = ^L me G log(m £ 2 //i 2 ) • (A7) 

7T 



In this limit the differential scattering cross section is given in eq. IS, while the full 1-loop 
result is obtained by multiplying it with \C^ /C^ L \ 2 . The numerical value for the ratio 
C^p /C^ L is shown in fig. [], obtained using the LoopTools package . 

The scalar-type DM-lepton interaction £ e fj = G(xT x x)(£cs£) with r x = (c§ + 20^75), 
induces DM-quark interaction through two-loop diagrams, see fig. 0. This contribution is 
relatively easy to compute in the limit of heavy leptons using the operator product expansion. 
First one integrates out the heavy leptons, thus matching onto the local dimension seven 
operator 

C eS = C s (ri — (xT x x)F^/e 2 , (A8) 

■mi 

where the Wilson coefficient is 

Cs = \<fic f s . (A9) 

This then enters a loop with two photons attached to the nucleon current. We evaluate this 
1-loop diagram in the heavy nucleon limit, which gives for the matrix element 

M = ^ — {u' x T x u x )Z 2 F(q) (u' N ±(l + l0 )u N ) , (A10) 



with k = \ / 2mNEd the recoil three- momentum of the nucleus, and F(q) the two-loop nuclear 
form factor. The resulting differential cross section is given in eq. |T9| . 

The derivation for scalar DM follows along the same lines. One matches onto the dimen- 
sion 6 operator 

C cS = C S ,M — (x ! X)F^/e 2 , (All) 



mi 

with the Wilson coefficient Cg^ — i^emCs^s, which then gives a matrix element for \N 
scattering 



C AB 2tt 2 k 



M = T7t4 —^Hq) + 7oM . (A12) 

lo7r mi 



APPENDIX B: DERIVATION OF THE COUNTING RATE FOR WIMP 
ELECTRON SCATTERING 

Using coordinate space Feynman rules we obtain for the matrix element for WIMP scat- 
tering on an electron bound in the atomic shell with quantum numbers nlm 

M r n t s ' = g J d 3 xVw(x) e*( k - k '- p ><rx«:r e x , (Bl) 

Here, k and k' are the initial and final momenta of the WIMP, and p' is the momentum 
of the electron in the final state. The momentum of the initial state bound electron is 
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determined by momentum conservation: p = k' + p' — k. Specializing now to the axial 
vector case T^ = T^ = j^j 5 we obtain for non-relativistic \ and electrons 8 

W^J 2 = \ E \MZ^\ 2 = ^mlmlG 2 \ij nlm (p)\ 2 , (B2) 

rr'ss' 

where the momentum space wave function ip n i m (p) is defined by 

</w(p) = f d 3 x Vnim(x) e- ipx = xm{p) Y lm {9, 0) , (B3) 



with the normalization 



(2tt) 3 

For the differential cross section we have 



<rp IUp)I 2 = i. (B4) 



dE d 



\M-nlm 


2 


32nE x E e v xe 


k + p| 



(B5) 



Here, E x ~ m x , E e w m e (using "C m e ), and t> xe is the relative velocity of the WIMP and 
the bound electron. The event rate is obtained by summing over all shells and integrating 
over the electron momenta in each shell: 

<K WES Po'fe f <Pp f M , , , „ daS s 



nlm 



dE d m x p dct J (27r 

The angular dependence of the wave function disappears due to the orthogonality relation 
Yi m (9, 0)Y"^(#, 0) = {21 + l)/47r for the spherical harmonics. In the laboratory frame 
the electron and DM momentum form an angle 9, so that cos 6 = kp/kp. For the integration 
over cos 9 from the d 3 p integral we have to take into account eq. which holds under the 
approximation E d <C m e < E e <C m x and implies 

cos9^- 1 -(^ + ^-) . (B7) 



v \ p 2m^ 

The leading order corrections to this expression give the kinematically available range for 
cos#: 



;min ~ 

m x vp 



2m e (E d - E B ) (m 2 y - 2m x E d ) . (B8) 



Then, |k + p| ^/m^v 2 — 2m x E d , and eq. |B6] leads to the expression for the counting rate 



given in eq. |30[ The approximations we have made in deriving the formulas are accurate up 
to about 10% for m x ~ 1 GeV, but the error decreases with increasing m x , and we estimate 
an accuracy of 0(1%) for m x > 100 GeV. 



Using the hydrogen-like atom approximation we esimate the relativistic corrections to be of order 20% for 
electrons from the Is shell of iodine and smaller for the other shells. 
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APPENDIX C: NUMERICAL EVALUATION OF ATOMIC MATRIX ELE- 
MENTS 



In this appendix, we describe how we compute the radial momentum space wave function 
x{p) for WIMP-electron inelastic scattering and the matrix elements {n'l'm'\e % ^~^ ^\nlm) 



defined in eq. 37 for WIMP-atom elastic and inelastic scattering. 



WIMP— electron scattering (WES). The momentum space radial wave function 
Xnlip), see eq. EH required for the evaluation of the event rate for WES, is obtained by 
splitting the coordinate space wave function ^ n ; m (x) into its angular part Yi m (9, </>) and its 
radial part R n i(r), and computing 



Xnl{p) 



An 



21 + 



:j-y^Vw(p) Yi m (6 p ,. 



4iri l drr 2 R n i(r)ji(pr) . 



(CI) 



Here, p is a momentum space vector with modulus p and arbitrary orientation (6 P , 4> p ), and 
Pi(cos8) is a Legendre polynomial. In the second line, we have used the orthogonality of the 
spherical harmonics, and in the third line, we have used Gegenbauer's formula |64| |, which 
relates the Fourier type integral over a Legendre polynomial to the spherical Bessel function 
of the same degree. 

The radial wave functions R n i(r) can be approximated by a linear combination of so-called 
Slater type orbitals (STOs) 



Rm(r) 



C-nifc " 



(2Z lk ) n ^ 



3/2 



y/(2n lk )\ 



(r/a ) nik - L exp(-Z lk r/a ) 



(C2) 



Here, Oq is the Bohr radius, and the parameters c n i k , ni k , and Z[ k are taken from |65 



With R n i(r) given in the form of eq. |C2| , we can evaluate |Cl] analytically, which gives 

,-i+n lk ( 27Ta o V /2 ( ipaoV (l + n lk + l)\ 



Xnl 



(p) = E- 



Cnik 2 



J lk 



x 2 Fi 



U2 + l + n lk ),l(3 + l + n lk ),l + l,- 



pa 



Hk 



(C3) 



with 2-F\(a, b, c, x) being a hypergeometric function. 



WIMP— atom elastic scattering (el- WAS). The matrix element for this case can be 
written as 



J2(nlm\e i( - k -^\nlm) = ^/ ^dOr 2 [^(r)] 2 ^(0, <j>) Y lm ( 



AKr cos 8 



(2/+1) ldrr*[RUr)} 2S -^ 



(C4) 
(C5) 
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with the abbreviation K = |K| = |k— k'|. This integral has the form of a Fourier Sine Trans- 
form, and can be evaluated efficiently using the Fast Fourier Transform (FFT) algorithm. 
For R n i(r), we use again the expansion eq. |C2| , with the coefficients taken from [J65| . 



WIMP— atom inelastic scattering (ie-WAS). Here the numerical evaluation of the 
atomic matrix elements is slightly more involved than for el- WAS because now Vvz'm'( x ) 7^ 
■0 n / m (x). We expand the factor e* Kx in eq. [IT] in spherical harmonics 0] and rewrite the 
angular integral over a product of three spherical harmonics in terms of the Wigner-3j 
symbols |J7|]. This gives 



(n'l'm'\e \nlm) = An drr 2 R nl (r) R n >i>(r 



Y,h{Kr)Y LM {QkAk) 

L,M 



V(2/ + l)(2/' + l)(2L 



1 


V 


L ) 


(I V 


L ) 


Vo 





o 


\m m! 


Ml 



(C6) 



where ]l denotes a spherical Bessel function of the first kind, and 9k, (Pk are the angular 
components of K. To compute the cross section, we need the expression 



I (n'l'm'l^lnlm) = {21 + 1)(2Z' + 1) J^(2L + 1) 



I V L 


2 



drr 2 R n i(r) Rn'v{r)j L {Kr) 



(C7) 



Here, we have used the symmetry and orthogonality relations of the Wigner-3j symbols 
and of the spherical harmonics. The expression in the second set of square brackets has 
the form of a spherical Bessel transform, which we evaluate by using an algorithm due to 
Sharafeddin et al. p5[ , based on rewriting the spherical Bessel function as a finite Fourier 
series, thus converting the integral to a sum of Fourier sine and Fourier cosine transforms. 
The initial state wave functions R n i are again given by eq. and ref. ||65|| , while for the final 



state wave functions R n 'i>, we use the hydrogen-like approximation (with an effective charge 
Z = 3 due to screening of the nuclear charge by inner electrons) since accurate tabulated 
wave functions for excited atoms were not available. Also, we consider only transitions from 
the Is, 2s, 3s levels to the first two unoccupied s shells; we have checked numerically that 
these transitions are the most important ones. Fig. |l^ shows some of the ie-WAS matrix 
elements for germanium and xenon. We find that for germanium (xenon) transitions from 
the Is (2s) orbital dominate. Our approximations should correctly reproduce the qualitative 
behavior of the matrix elements at large momentum transfer, and should lead to a good order 
of magnitude estimate of the ie-WAS event rate. 



APPENDIX D: ON THE EQUILIBRIUM OF DM CAPTURE AND ANNIHILA- 
TIONS IN THE SUN 

In the calculation of the neutrino flux from DM annihilations in the Sun we have assumed 
that WIMP captures and annihilations are in equilibrium, which makes the result indepen- 
dent of the DM annihilation cross section (a ann v). Here we comment on the validity of this 
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FIG. 10: A few examples for ie-WAS matrix elements of germanium and xenon. The plot shows 
transitions from the Is resp. 2s orbitals to some of the lowest unoccupied states. 

assumption. Let us first estimate the cut-off scale A for the effective theory description of 
the DM-lepton coupling. For the two examples of V <8> V and A® A couplings, the neutrino 
bounds are of order <r° e ~ 10~ 43 cm 2 and 10~ 38 cm 2 , respectively, see figs. ^ and [7[ From 



eq. [13] we can estimate the corresponding cut-off scales as Ay ~ 100 GeV and A^ ~ 10 GeV, 
where we took coupling constants cf to be of order 0(1). In DM annihilations the four- 
momentum transfer squared is of order m 2 . For m x ~ 10 GeV, relevant for WNS, the 
WIMP annihilations may then also be described by effective field theory. Using effective 
interactions in eq. [I] (extending them to neutrinos), we find 



G 2 m 2 „ m 2 „, „ , / a 



Vector: (* m v) ~ = ^ ~ 10" 24 cm 3 s" 1 ( U ^ ) . (Dl) 



7r m 



In the WES case, however, the effective theory typically cannot be applied since the mo- 
mentum transfer for annihilations is above the cut-off scale. Therefore, in general we cannot 
make model independent statements about (cr ann v) without specifying the UV completion 
of the effective vertex. An order of magnitude estimate can still be obtained from dimen- 
sional analysis as 

4 

9 1 rv-21 ^3 „-l v „A f m X 



Axial: (a ann v) ~ ^ ~ 10^ cm 3 s^x/ — - S— , (D2) 
mi V 100 GeV/ 

A, 

with g a typical coupling constant between leptons and the dark sector. 

Equilibrium of WIMP capture and annihilations is obtained if tanh 2 (t /r) is close to 1, 
see eq. [|3| Fig. |TI| shows the values of (cr ann t>) for which t Q /r = 1 and 5 as a function of 
m x . The values of scattering cross sections a xe for V <S> V and A® A Lorentz structures 
were chosen to be above (but close to) the Super-Kamiokande bounds shown in figs. [| and 
|7|. Since tanh 2 x ~ 1 for x > 5, WIMP capture and annihilations are in equilibrium in the 
Sun for values of (u^v) above the curve for t /r = 5. Comparing eqs. [Dl] and [DD with 
the ranges shown in the figure we conclude that the assumption of equilibrium is very well 
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m x [GtW] 

FIG. 11: Contours of io/r = 5 and i©/r = 1, see eq. For the case of vector (axial vector) coupling 
we have used a scattering cross section of cr° e = G 2 ml/n = 10~ 43 (10 -38 ) cm 2 , motivated by the results 
of the Super-Kamiokande bound. For values of (a ann v) above the curve for £q/t = 5, WIMP capture and 
annihilations are in equilibrium in the Sun. 



justified in the cases of our interest. 
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